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We investigate dynamics of large scale and slow deformations of layered structures. Starting from 
the respective model equations for a non-conserved system, a conserved system and a binary fluid, we 
derive the interface equations which are a coupled set of equations for deformations of the boundaries 
of each domain. A further reduction of the degrees of freedom is possible for a non-conserved 
system such that internal motion of each domain is adiabatically eliminated. The resulting equation 
of motion contains only the displacement of the center of gravity of domains, which is equivalent 
to the phase variable of a periodic structure. Thus our formulation automatically includes the 
phase dynamics of layered structures. In a conserved system and a binary fluid, however, the 
internal motion of domains turns out to be a slow variable in the long wavelength limit because of 
concentration conservation. Therefore a reduced description only involving the phase variable is not 
generally justified. 

PACS number(s): 68.10.-m, 68.60.Dv, 68.15,+c 
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I. INTRODUCTION 



Layered structures appear in various systems in nature. Smectic liquid crystals ffl, mcsophases in a microphase 
separation of block copolymers S, magnetic bubbles in magnetic thin films || and condensed phases of Langmuir 
monolayers , || are typical examples in thermal equilibrium. 

A layered structure can also be observed in pattern formation far from equilibrium Q. Turing was the first to 
predict a stable periodic structure in a reaction-diffusion system. Recently, a Turing instability has been observed 
experimentally in a chemical reaction in a gel [Q. Other well known examples are roll patterns in Rayleigh-Benard 
convection and Williams domains in the electrohydrodynamic convection of liquid crystals. 

In order to study dynamics of these structures, a coarse-grained description is often employed. Suppose that there 
is a periodic layer with period £ in x-dircction. This is expressed by a periodic function X(x + <jy) = X(x + £ + (f>) 
where X is a physical property and is a phase and is arbitrary in a uniform extended system. When the periodic 
structure is slightly deformed, its large scale slow motion can be described by allowing to depend on x, y, and t. 
Since is a neutral mode, it is one of the slowest variables in the long wavelength limit. The theory based on this 
idea, which is called phase dynamics, was developed in pattern dynamics far from equilibrium §,0. 

An essentially similar method can be applied to layered structures in equilibrium jllj. One can express the free 
energy increase due to deformation in terms of the local displacement field as in continuum elastic theory. The local 
displacement field corresponds to the phase variable and the dynamic equation for it can be written down by using 
the free energy functional. 

Many of the theories assume a sinusoidal function (and only a few higher harmonics) of X(x) in the derivation of the 
phase equation. This is valid when the system is in the vicinity of a bifurcation (critical) point where a uniform state 
loses stability and a periodic structure appears supercritically. However, when the function X changes abruptly at the 
interface separating two adjacent layers, we need to develop an alternative method for dynamics. Phase dynamics is 
still meaningful in this case for long wavelength deformations, but the derivation of the phase equation is not possible 
using available theories. 

In this paper, we address the above problem. Some years ago, Kawasaki and one of the authors 12 derived a 
free energy functional of general periodic structures in terms of the displacement fields without assumingthe slow 
variation of X. Recently a similar analysis has been applied to a model specific to Langmuir monolayers Jjjfl . Here, 
on the other hand, we are concerned with dynamics. The static aspects can be obtained as a byproduct by taking 
the equilibrium limit. It is emphasized that our theory does not necessarily rely on the existence of a free energy 
functional or a Lyapunov functional. As a consequence, we can deal with the systems both near equilibrium and far 
from equilibrium in a unified way. 

Our strategy is to derive the interface equation of motion for each domain boundary. Suppose that X varies within 
one period as X = 1 for < x < 2w and X = — 1 for 2w < x < I where w(< 1/2) is a constant. Thus there are 
interfaces at x = and x = 2w bounding a domain where X — 1. We assume the interface width is infinitesimal 
compared to w and £. It is generally expected that these two interfaces are coupled to each other as well as to those 
in other domains. We can derive the phase equation of motion from the interface equation by retaining only the 
average displacement of two interfaces at each domain. The other degree of freedom which makes the center gravity 
of a domain time-independent represents the internal motion of domains. In some situations, internal deformations 
of each domain are relevant to the stability. The present theory covers such a domain instability. 

The phase variable describes the Goldstone mode associated with the translational symmetry breaking due to a 
periodic structure in an extended uniform system, and hence it is a slow variable in the long wavelength limit. In the 
terminology of lattice vibrations this corresponds to an accoustic mode. On the other hand, the internal motion of 
domains is similar to an optical mode and has a finite energy in the limit of long wavelength deformation. As will be 
shown below, this gives rise to a finite relaxation rate of the internal mode for small wavenumbers in a non-conserved 
case. Therefore one may eliminate adiabatically the internal mode. However, as we will show, the situation is entirely 
different in a conserved case and a binary fluid. When the interface width is infinitesimal, and there is no mass flux 
through the system boundary, concentration conservation is equivalent to a conservation of total area (volume in three 
dimensions) of domains. Therefore, the internal mode should also be a slow variable in the long wavelength limit. 
Our main concern in this paper is to clarify how the phase dynamics are affected by the extra slow mode originating 
from conservation. 

We start with a time-dependent Ginzburg-Landau (TDGL) type equation describing the system at the coarse- 
grained level. A general method is available to derive the interface equation of motion jF4| from the TDGL equation. 
In binary alloys or in crystal growth, on the other hand, it is also well-known that the same equation can be obtained 
by starting at the macroscopic level with the simple diffusion equation for the minority atoms and employing the 
Gibbs-Thomson relation as the boundary condition at an infinitely sharp interface |lq] . Thus one may wonder why 
the TDGL approach is necessary. There are, of course, several reasons. First of all, most simulations of spinodal 
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decomposition and other phase separation phenomena in binary systems use TDGL equations since for many purposes 
it is more efficient than the macroscopic approach. Hence the interfacial analysis based on the TDGL equation is 
necessary for a complete understanding of the simulation results. Secondly, the interfacial method is valid beyond the 
simple diffusion equation which has been used only for macrophase separation and for crystal growth. As mentioned 
above, we will show that our method can be useful for both thermal equilibrium structures and dissipative structures 
out of equilibrium. 

The organization of this paper is as follows. In section 2, we make a general argument on layered structures in 
equilibrium. In section 3, we derive the interface equation of motion for a periodic structure in a non-conserved case. 
In section 4, we verify that the vanishing of one of the diffusion constants in the phase equation is equivalent to free 
energy minimization with respect to the spatial period. The method is extended in section 5 to a conserved system. 
In section 6, we treat a layered structure in binary fluids. The results obtained are discussed and summarized in 
section 7. Some of the details of the method used in section 3 are described in the Appendix. 



II. PERIODIC STRUCTURES 



A typical example of a layered structure in a system with a non-conserved order parameter is a magnetic thin film 
p6[ . Experiments fL7| show stripe and bubble structures of magnetic domains depending on the magnitude of the 
external magnetic field normal to the film. The free energy of this system can be written in terms of the local order 
parameter it(r) as 



where 



and 



F = F GL + F LR (1) 



Fgl = I dr\-{Vu) 2 + W(u) - hu] (2) 



F lh = \\ drj d?G(r, ?)[u(f) - c] W) - c] (3) 

Note that the free energy is generally divided into two parts. Fgl is the short ranged part and is the usual Ginzburg- 
Landau form with W(u) a double well potential in the ordered state and h the external field. The parameter e, which 
is a measure of an interfacial width, is assumed to be sufficiently small. The other part Flr contains the long range 
interaction, where a is the interaction strength and c is a constant. The form of the kernel G(r, r 1 ) depends on the 
system considered. For instance, it is given in a thin magnetic film of a strongly uniaxial material by Jt6| 

f Air 

G{r,r') = J dq—[l-exp{-qD)]exp(iq ■{?-?)) (4) 

where the integral over q is taken in a two-dimensional space parallel to the film surface. Note that the order parameter 
u in this case is the magnetization normal to the film which has thickness D. The spatial period of a modulated 
structure is larger than D, and the interface width e is much smaller than D, i.e., e << D << I. Probably the 
simplest choice of G(r, f*) which provides us with a periodic structure in equilibrium is the Coulomb Green function 

-V 2 G c {r,r') = 5{f-r') (5) 

This type of long range interaction appears in the theory of phase separation of block copolymers |l8| . The constant 
c in (|3j) is equated to the spatial average of u to avoid a divergence in the long range interaction. 

In terms of Fourier components, the sum of the short range and the long range interactions contains the quadratic 
part 



1 

2 

where 



r(g)w<fu-(j (6) 



T{q) = e 2 q 2 + aG q - r (7) 
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with G q the Fourier transform of G(f — r"). The constant r comes from the bilinear term in the local part W(u). 

In order for a stable periodic structure to exist, we assume that T has a minimum for a finite value of q and that 
the free energy (|l|) takes the lowest value for such a modulation. Hereafter, to be specific, we consider the case where 
G(r, r") is given by the screened Coulomb interaction 

— (V 2 — K 2 )G(r, f 1 ) = 5(r — f) (8) 

where 1/k is the screening length which is assumed to be much larger than e. In this sense, the interaction G(r, r 1 ) 
is still long ranged. The interfacial theory for other forms of G(r, r 1 ) can be constructed similarly; in particular, the 
pure Coulomb case follows by taking the n — > limit of the interface equations. 



III. NON-CONSERVED SYSTEM 

The time-evolution equation of a non-conserved order parameter u is determined by the free energy functional (jl]) 
according to 

*t--*L (9) 
dt Su 1 ' 

It is noted that the right hand side contains the nonlocal long range interaction. When we derive the interface equation 
of motion for layers, it is more convenient to express it in a local form. To this end, we write (^|) in the form of a 
coupled set of equations as 

Tl ^= e ^u + f(u)-v + h (10) 



dv(r) 9 

t 2 ^- l = V v - K 2 v + au (11) 
at 



where we have introduced two parameters t\ and ti and /(u), defined by 



/(«) = (12) 
au 

The variable v in ( [Til) is an auxilliary field and has no physical meaning for the model (|9|). However the set of 
equations (10) and ( |il|) appears in pattern formation far from equilibrium as will be mentioned below. Without loss 
of generality, we may put the diffusion constant for v equal to unity. In the limit t-2 — » 0, cq.(p"T|) can be solved for v 
by using G(r, f") defined by (||). Substituting it into eq. ( [To| ) we recover eq. (||). 

For a general form of the kernel G(r, r 1 ) in the long range interaction, we may replace ~S/ 2 v — k 2 v by J dr f T(r, r , )v(f f ) 
where r(r, r") is defined through the relation 



dr"G(r,r")T(r",r') = 5{r-r") (13) 

It is remarked that eqs. (|l^) and ([ll]) are model equations for pattern dynamics far from equilibrium. For instance, 
the Belousov-Zhabotinski reaction in the excitable regime can be described, in its simplest approximation, by eqs. 
( |l0| ) and ([ll]). In fact this set of equations has been extensively studied during the past two decades by applied 
mathematicians p^ , |2(J. It is well known that eqs. ( |Toj ) and (jll]) undergo a Turing instability when e is small so 
that a periodic structure can be formed. To our knowledge, however, an interfacial approach has not been applied to 
a periodic solution except for ref . |^lj . 

We solve eqs. (|l(]) and ( |IT|) in the limit e — > to derive the interface equation of motion. We assume that when the 
diffusion terms are absent, these equations have two stable uniform solutions (u\,vi) and (u2,V2)- As shown in Fig. 
1, we consider a layer of domains in two dimensions where the regions u = u\ and u = 112 are arrayed alternatively. 
We assume the width of a domain in which u = u\ is equal to 2w and the period of the layer is I. 

Now suppose these striped domains are gently curved as in Fig. 1. We introduce the deformations CniUy ^) (Cn (2/> £)) 
around the equilibrium position of the right (left) interface in the n-th domain in which u = u\. Our aim is to derive 
the linearized equation of motion for . For concreteness, we hereafter assume the following form of f(u) 

f{u)= l -u{l-u 2 ) (14) 
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As was mentioned above, the parameter e is a measure of an interface width. Suppose that we are concerned with 
the motion of a specific interface. In the limit e — > 0, the density u varies rapidly through the interface, but the spatial 
variation of v is much slower. Therefore we may replace v in eq. ( |To| ) by its value at the interface, which is denoted 
by vi and which may depend on time and position on the interface. Thus eq. ( [To|) can be written as 

T 1 ^ = e 2 V 2 u + f(u)-v I + h (15) 

This is just the TDGL equation with the "external field" , vi — h, and hence the interface equation of motion for this 
system is readily obtained: 

t\V = e 2 H - 3e(vi - h) (16) 

where, for simplicity, we have assumed that the absolute value of i>j — h is sufficiently small so that the gap of the 
density u at the interface is approximately equal to 2. V is the normal component of the interface velocity and H 
the mean curvature. These are defined as follows. First let us introduce a field variable 4>(f,t) which is positive 
(negative) in the region u > (it < 0) and ip(r,t) = specifies the interface configuration. The velocity V and the 
mean curvature H are expressed respectively as 

|VV>| dt K ' 



H = V-n (18) 

with 



(19) 



Note that n is the inward normal of a u > domain. 

From the above definition of ip, deformations may be written as ip — —x + ( + (y,t) and ip = x — £~(y,t). 
Substituting these into ([17]), (18) and (|l6|), we obtain up to linear order in 



TLdfi/dt - e 2 8 2 Ct/dy 2 T 3e( V/ - h) (20) 
The last term indicates that a domain where u > shrinks as Uj — h is increased. 



The remaining problem is to determine vi as a functional of so that eq. (20) becomes a closed set of equations. 
In what follows, we need to solve eq. ( p"l| ) for a given set of C^- This was carried out in a previous paper [|T) for 
general values of T2 to study a dynamical instability of layered domains. Here we consider only the limit T2 — > 0. In 
this limit, the solution is given by 

v(r,t) =a J dr'G(r,r l )u(r / ,t) (21) 

When the absolute value of Vj — h is sufficiently small as we have assumed, the uniform solutions U\ and u-i may be 
approximated as u\ = 1 and ui = — 1. Actually only the difference u\ — enters in the theory given below. The 
correction is however, of order (vi — h) 2 , i.e., u\ — U2 = 2 + 0((vj — h) 2 ) and hence does not affect the first order 
correction in (|l6|). Therefore eq. ( ^l| ) can be written as 

dx'^-^expiiq- (f- ?))] (22) 



where we have defined 



tfa=ml±w + £ (23) 



5 



The value Vj in the n-th domain, which generally takes a different value at the right and the left interfaces, is given 
by 

vf = v(b± >y ,t) (24) 

Now we calculate vj to first order in £ . Some further details are provided in the Appendix. The zeroth order 
solution is given from eq. ( p3| ) in the Appendix by 

where reciprocal lattice vectors are given by Q = 2im/l with n integers (n = 0, ±1, ±2, ...). The first order solution 
is given from eq. ( |94] ) in the Appendix by 

= -2aC(y) J £ Q^T^t 1 cos2 ^) 

Q 

+ °J 2ir e J 2n£^ (k-Q)*+p 2 + K 2 

xlti p -e- 2iQw tkJ (26) 



9tt I Ok P Z—i 



2tt 7 2tt ^ ^ (fc - Q) 2 + p 2 + k 2 
x [e^C^-Q,] (27) 



In the derivation of eqs. (26) and (|27|) we have introduced the Fourier transform 

f dk f dv 

£(y)=J^J etPMnt ±*>) + iyp] (28) 



and used the Poisson formula 

oo 1 

£ exp(ifcnO = -£*(fc-Q) (29) 

n=— oo Q 

Substituting these into eq. ( p0[ ) we obtain the interface equation of motion. The zero-th order solution is particularly 
important, which reads 

^ = 7?^-^ (30) 

This specifies the ratio w/£ of these unknown parameters. It is readily verified that when h — 0, w/i? = 1/4, as it 
should be for the symmetric case. Equation ( p0[ ) is equivalent to dC^/dt = without deformations. Hence this is an 
equilibrium condition for a nonconserved system. 

From the first order equation, we obtain the coupled set of linear equations for : 

Q 

+ — £ (fc-Q) 2 + P 2 + « 2 + e_2iQ " c fe,p] ( 31 ) 



G 



d 2 2 6ea ^-^ 1 

n ^ = " [ep " — Z.q^t^ (1 " cos2Qw)]c ^ 

Q 

+ — E (fc - Q)2 + p2 + K 2 i e2iQW C P ~ Cfc, p ] (32) 

Q 

In order to analyze the motion of long wavelength deformations, we make the following transformation 

ti ] (y) = l(ti(y) + Q(y)) (33) 
C ) (y) = l(C(y)-Q(y)) (34) 

In real space, (n^ describes the deviation of the center of gravity of each domain while (jff 1 is the internal deformation 
of a domain. Note that (ffl ^ but Q, = for a uniform translation of a domain. Thus (n^ is a phase variable. 

Since we are concerned with large scale deformations of a modulated structure, we may apply a gradient expansion 
to eqs. (|3l| ) and (^). That is, we expand the coefficients of these equations in terms of k and p. Retaining the lowest 
nontrivial order, we obtain from eqs. (J3l|) and (p2[) 



|c2 = -W + !=<«■ - i) Y. why-' 1 - » s2<3 "*'S 



12iea ( 2 ) v-^ Q 



e - fc €E (Q2 + K2)2 sin2 ^ 



4« = -^ + TS^?-^KS 

Q 



Note that the right hand side of eq. (35) vanishes in the limit k,p — > whereas that of ( |3q ) remains finite in this 
limit. Hence we may assume relaxes more rapidly, and we may put dC^ /dt = in d36|). We may then substitute 
the solution into eq. ( |35"| ) to obtain 

T 4 C S = ~x [i ^ 2 + D ^ 2 + Kpi]c % (37) 



where 



D - = 4 i- 8 iJ2whr {1 - cos2Qw) (38) 

Q 

24a 1 , , 

D W = — L ( q 2 + k2)2 ( 1 - cos2 ^) 

16a l^Qsin2Qw 2 



and 
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K = ^E w h^ il - coa2Qw) (40) 

In eq. ( p7h, w e have factored out 3e£/A. The reason will become clear in eq. ( |50| ) below. Eq. (|37]) with the coefficients 
given by (|38;), ( j39[ ) and (flO]) is the general form for layered structures. The coefficient K is always positive and Dy 
is positive at least when w w £/4. Thus the layered structure in the present model system does not cause any linear 
instabilities in this regime. 

The coefficient D±_ must be identically zero when the spatial period t is the equilibrium period. This will be verified 
separately in the next section. 



IV. DETERMINATION OF THE PERIOD 



When the system has a free energy functional, the equilibrium period should be determined by its minimization. 
In this section we show that this is identical to the condition D± = 0. Although this is a well-known fact (see, e.g., 
ref. Jl^l ), a verification is necessary to confirm the validity of the present approach. 

We consider a one-dimensional system with the system size L(» €) and write the free energy per unit length as 

F = i I dx [ dx's{x - x')u{x)u{x') + ~ / dxW(u) (41) 



2£ . 

where s(x — x') is the nonlocal interaction. In the case of ([!]) it is given by 

B 2 

s(x-x') = e z -^^8{x-x') + aG(x-x') (42) 

The second term of ( ffll) is the local part. We assume that the solution of 5F/5u — is spatially periodic. 

Now suppose that the system is expanded uniformly so that the period becomes 1(1 + 8) (S << 1). Thus u and L 
are transformed as u(x) — * u(x/(l + 8)), L — > L(l + 8). The free energy changes as 

F = ^Jdyj dy's((y - y')(l + S))u(y)u(y') + i J dyW(u) (43) 

where x/(l + 8) — y. The free energy increase AF is, therefore, given up to 0(8) by 

AF = ^- J dy J dy' s(y - y')u{y)u{y') 

^J^jmv-y')[^^Hy)<y') 

= -^lj d yj dy'{y-y')s{y-y')[^]u{y') 

This may be rewritten in terms of the Fourier components as 

If u is the equilibrium solution, AF must be identically zero. 

When s(x — x') is given by ([42|), s q = e 2 q 2 + a/(q 2 + k 2 ) so that we have 

Let us calculate (^6|) explicitly. First of all, we define the equilibrium interfacial energy a. When f(u) is given by 
(|l2|), it is calculated as 

(47) 



(44) 
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where the integral domain over x is restricted to the vicinity of an interface. Note that the interface width dependence 
of a is different from that in ordinary critical phenomena. However, this is only apparent and is due to the definition 
of the free e nerg y (||) where e appears as a coefficient of the gradient term. 
By using (^) we n °te 

e 2 f dq 2 e 2 f du. 2 cr 



L J t^ U « U -« = L J ^ = 2 7 ^ 

The factor of 2 comes from the fact that there are two interfaces in each domain. 

Since we have assumed that u = 1 for < x < 2w and u = — 1 for 2w < x < £ away from the interfaces, we obtain 

a f dq q 2 . 8dv-^ 1 , . , s 

l J 2^WT^Y )UciU - (1 = T ^ WT^y [l - cos2Qw) (49) 

Putting this together with (|48| ) we finally obtain 



This precisely agrees with D± — given by (38) 



V. CONSERVED SYSTEM 

In a phase separation of binary mixtures, the order parameter is a local concentration and is a conserved quantity. 
The time-evolution of a conserved system is modelled by 

This equation with the free energy given by (0) was introduced to study a microphase separation of block copolymers 

& ' _ ' 

When the long range interaction is absent in (||) , eq. ( pl| ) is called the Cahn-Hilliard equation for spinodal decom- 
position, and the interface equation of motion has been derived. However, an interfacial anaylsis for multi-layered 
systems given by (^lj) has not previously been obtained. 

First we reformulate the diffusive model ( pi] ) by introducing an auxiliary field v. 

Ti Mn =e 2 v 2 u + f,j_ a [ df ?G(r,T<)u(r*,t)-v + h (52) 
at J 

9v(f) _ 2 . du(r) 



where the long range interaction has been put in ( p2[ ) . If the constant a were zero in ( |52|), the above set of equations 
would be essentially the same as the so-called phase field model for melt growth p3|, |24j]. In this case, the variable 

v is the local temperature and the last term in ( |53| ) stands for the latent heat production upon crystalization. 

Eq. d5l| ) can be recovered as follows. We operate V 2 on both side of eq. (|52]). Putting t 2 = in eq. (p^), we may 
write V 7 ^ in terms of u, and hence (|52|) becomes a closed equation for u. Taking the limit T\ — > 0, we obtain eq. ([5l]). 
Thus the expression of the diffusive model as eqs. (|5^) and ( p3[ ) enables us to analyze the interface dynamics almost 
parallel to the non-conserved case. 

However, there are two important differences. One is the fact that Uj now contains the time derivative of because 
of du/dt in (|53|). Putting n = 0, eq. ( |f5| ) becomes 

= e 2 H - 3e(uj - h) (54) 

It should be noted, however, that the contribution from the third term in eq. ( |52"| ) is also included in vi. That is, 

vi consists of two parts vi — v^p + v^p where comes from the third term in ( |52| ) and is given by (|2^) and j27| ) 



whereas w^ 2 ' comes from the fourth term in ([52|), which will be evaluated below. Eq. ( |53| ) can be solved as 
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v(f,t) = J d fG e (f,f)^^ (55) 

It should be noted that the Green function in ( p5| ) is not the screened one (||) but the Coulomb Green function (||). 
Using the fact that 

^ = -2j2lLd(x -m£ + w- C"(y)) - C+S(x -ml-w- C+(y))] (56) 
m 

where the factor 2 comes from the gap of u at the interfaces, we readily obtain up to linear order in 

where ?)j is the Fourier component of v\ . In these expressions we have emphasized that the n — component 
is included in the summations. 

The other difference from the non-conserved case is that we have to take account of the conservation law in the 
interface equation of motion. In the non-conserved case, we consider the zeroth-solution separately as (|30| ) which 
determines the ratio w/£. However, in a conserved system this ratio is fixed uniquely in the strong segregation limit 
by the average volume fraction. In fact, when there is no current across the system boundary, the conservation law 
requires 



d 
dt 



]T J dyt^ = o (59) 



The constant h in eq. (|52) plays a role of a Lagrange multiplier associated with the condition (|5£ 
Substituting ( |57| ) and (pq) into (^4j) and retaining the zero-th order term, we obtain 

7 E (k - Q) 2 +p 2 [ (l-cos2Q W )Cg-isin2Q<g] 



Q(=o) 



7 E {k-Qy+f [ ( 1 + cos2 ^)CS + l sin2Q<S] 

r 9 n 12ea x— \ 1 _ . jo\ 

= -[eV + — £ cos2QHCS 

Q 

- lg F fc cSE (Q2 + K2)2 ^ 2 ^ 

Q 

- CS(k)S( P ) (61) 

where 
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2a x - sin 2Qw 
T Q(Q 2 + k 2 ) ~ 



(62) 



The right hand sides of ( |6C| ) and j6l| ) have been expanded in powers of k and p and are the same as those in (|35| ) and 
( |36| ) respectively. The unknown constant h or, equivalently C, is determined by (59). That is, in the limit k,p — > 
in eq. ( |6l| ) we obtain 



C 



12ea v cos 2Qw 



*■>(), 



(63) 



Note from (|63|) that the "zero-th" order term is actually not strictly zero-th order, but it turns out to be linearly 
dependent on C} 2 ^ . 

It should be noted that the summations on the left hand side of eqs. ( |60| ) and (61) contain the Q = component. 
In particular, eq. ([n]) produces the factor l/(k 2 +p 2 ) which originates from the conservation law and is dominant in 
the long wavelength limit. Thus eqs. (pOh and (O) are approximated as 



Sfc,p 



-r« l(D ±P 2 + Dk 2 + Kp*)(£l - ikB 2 cfl\ 



(64) 



eS = -rW(* 9 +p a )[fl 1 cS+i*fl ! icS] (65) 

where 

24a ^ 1 - cos 2Qw 

D = —\ {Q 2 + n 2 ) 2 (66) 



r(1) = \ (67) 

o 1 — cos 2Qw x ' 

2^Q=£0 Q2 

T (2) = y Q (68) 

Bl = ?y^ (69) 
Q 

16a ^ Qsin2Qw 

52 -— Z. (Q 2 + k2) 2 ( 7 °) 

We have omitted the C term, and D^and K are given by eqs. ( |38[ ) and (^), respectively. The factor has been 
factored out such that the term in (^4|) and ( |65| ) have the same coefficient. In this way we have a coupled set of 
equations for Q l ) and ^ 2 \ whose relaxation rates vanish in the limit k,p — * 0. 

There are two cases where the relaxation of C*- 2 -* is much faster than that of £W. One is the limit n — » where Si 
becomes very large because of the Q = component. Note that D remains finite. The other is the case that D± = 
and k 2 « p 2 . In these situations we can eliminate <^ 2 ) by putting Q 2 ' = 0: 



|<£ = -T^[D\\k 2 + Kp*]$> (71) 



where Dii is the same as in (£ 
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VI. BINARY FLUIDS 



Equation (|5^) has been applied to a microphase separation of diblock copolymers. When the sizes of the two blocks 
are nearly equal, the system exhibits a lamellar structure in the microphase separated state. The kinetics of this 
phase separation have been studied mainly by omitting any hydrodynamic effects. However, in some situations such 
as when the system is subjected to shear flow, hydrodynamic effects are known to cause an instability of a lamellar 
structure (25). 

Another example in fluids is the formation of a periodic structure within a Langmuir monolayer P6| , [ gfl . However, 
since this is an adsorbed system between an air/water interface, the hydrodynamics are different from that in bulk 
p8| . One needs a separate treatment of its interface dynamics which we do not enter here. 

In this section, we consider the conserved system coupled with the local velocity field v. 

^-+v- Vu = V 2 he 2 V 2 U -/( U ) + 0] (72) 

fj?) — * 

— = V V 2 v + ( Vu)[-e 2 V 2 u -f{u)+<t>] (73) 

V 2 0- K 2 <f) + au = (74) 

where rj is the shear viscosity. In order to avoid confusion of the velocity field v with the auxiliary field v, we have 
introduced the notation <p instead of v for the auxiliary density. We have imposed the incompressibility condition 
V • v = and eliminated the pressure term. Thus v has only transverse components. The nonlinear convective term 
v ■ Vv has been ignored in eq. (|73| ) . 

A general form of the interface equation can be derived as follows. We temporarily ignore the diffusive coupling 
of the order parameter given by the right hand side of ([72]) and consider only the hydrodynamic coupling v ■ Vu. At 
the final stage we will take account of the diffusive effect given by the right hand side of eq. ([72]). Thus the interface 
velocity V is given from ([72]) by 

V = -vi ■ ft (75) 

where the unit normal vector ft is directed inward to a u = 1 domain and vj is the fluid velocity at the interface. 
Putting dv/dt — 0, as usual assuming velocity fields relax rapidly compared to interface motion, we can solve eq. 
@. Substituting it into @, we obtain 

V(o) = J da'^n a {a)T a P{a,a')nP{a')[aH{a')-2a(t>i{a')] (76) 

where a and a 1 stand for the position on the interface. The tensor T a " is the Oseen tensor whose Fourier transform 
is given by 

Tf = \[^- q -^-\ (77) 

The constant a is the interfacial energy defined by ([47]). The value is the same as vi given by ( p6| ) and (p7|). In 
this case the zero-th order constant part of (f>i does not contribute to (|7^) because of the relation 

f da' n ° (a) TQ/3 {a 1 a')nP{a')=0 (78) 

a, /3 

Hence we may consider only the first order corrections given by (^6|) and ( p7| ) for 4>i in ([76]). 

The interface equation of motion can be derived from (f76|). Note that ft has only the x-component to lowest 
order and n x — — 1 for a right interface and n x = 1 for a left interface. Using the facts that V — idC^/dt and 
H = ±d 2 ( ± /dy 2 , we obtain 

Q(=0) Q' ^ 
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where £ is a column vector with £1 = and (2 = Cfep an d M{Q) ^ s a 2x2 matrix with the components Afn = 
1 — cos2Qw, A/12 = — is'm2Qw, A/21 = isin2<3u> and M22 = l + cos2Qw. The 2x2 matrix N(Q) has the components 
Nu = 1 + cos2Qw, A^i2 = zsin2Qw, A^2i = —isixi2Qw and N22 = 1 — cos2Qw. 

Now we expand ( |79| ) in powers of k and p where it should be noted that the summation over Q contains the Q = 
component. Using the relations 

k ' p v (k 2 + P 2 ) 2 y ' 

p 2 

T k-Q, P « ^4 (8!) 



eq. ( |79| ) becomes in the long wavelength limit 



r\(k 2 + p 



CS = -^[5iCg+i*B 2 Cg] (83) 



where 



1 v-^ 1 — cos2Qu; 
^=2^ nT^ ( 84 ) 



The reason for the appearance of the factor p 2 in (83) is clear. We have so far ignored the diffusive effect in ( |72|) and 
considered only the hydrodynamic interactions. Therefore, the concentration must be conserved within each domain, 
and hence the effective diffusion exists only along the y-axis in (|S3"|), i.e., along the interface. 

By comparing ( |82"| ) and (|83|), we note that if k 2 << p 2 and if << jBi for = 0, the relaxation of C' 2 ' is rapid 
so that we can eliminate it to obtain 

s^-as^W+'WKS (85 > 

This is also justified when k — > since the constant Bi becomes infinite in this limit. The coefficients are the same 
as those in the previous sections. If necessary, we may add to (pa) the contribution from the diffusive effect given by 
eq. ([n]). This equation has been obtained by several authors [ |29| , fjCfl . However, under general conditions, the phase 
variable has to be coupled with the other slow mode as in the set eq. ( |32"| ) with eq. 

VII. DISCUSSION 

We have presented a systematic method to develop the equations governing interface dynamics in layered structures. 
It is emphasized that the basic equations (|l0|), ( |TT| ) and ([52]), ([53]) do not possess any Lyapunov functional when t\ 
and T2 are finite. Although we have taken a limit in which these parameters vanish, this is simply because we have 
demonstrated our theory mainly for the systems near equilibrium. The interface dynamics can be formulated without 
taking these limits as has already been shown for a nonconserved reaction-diffusion system in Ref . |2~]}] . In the present 
paper, we have extended the methods to conserved systems and binary fluids. 

The motion of domains can be divided into two parts. One is the collective motion of center of gravity of each 
domain, which is described by a phase variable . The other is the internal motion of each domain denoted by . 

Equations of motion for £W and £*- 2 ) in conserved and non-conserved systems take the following form: 

*^ = -r«-** (86) 
at <5C (1) 
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where a = for a non-conserved case and a — 2 for a conserved case and H (which should not be confused with the 
mean curvature defined by (|l8|)) is given by 



H=\J d X dy[D( d -j-f + D^f + if (^) 2 + B l( C (2) ) 2 + 2S 2 ^ C (2) ] (88) 

This Hamiltonian is invariant under the transformation x — > — x since C (1) chan ges its sign whereas C*- 2 "* is invariant. 
Using the same reasoning, the coupling term (dCp^ /dy)^ 2 ^ should not exist. In the case of binary fluids, we can write 
the equations formally as 

8^ _ 1 8 2 SH 



dt 7y(V 2 ) 2 dy 2 5(W 
dt j] dy 2 5(( 2S > 



(90) 



These sets of equations show that although has a "mass" term with coefficient B\ in H, it does not produce a 
finite relaxation rate of C*- 2 -* in the long wavelength limit except in the non-conserved case. In the conserved system 
and binary fluids, the variable is generally a slow mode in the limit of long wavelength deformations because of 
the conservation law. As we have remarked, £( 2 ) can be eliminated only when the deformation normal to the layers 
is spatially sufficiently weak compared with the undulation of the layers. In fact, eqs. (^9|) and ( p0| ) give us the 
relaxation rate fr 2 ' of ^ 2 ^ approximately as 

(2) _ 7 ^ x Bl k 2 

V 1 jB 2 (k 2 +p 2 ) 21 [ ' 

where the second term is the correction due to the coupling with which is, however, not necessarily small for 
p -> 0. 

Note that the case of pure Coulomb long range interactions has to be handled carefully with the limit k — > 0. In all 
cases studied here, when k is sufficiently small, the internal mode Q 2 ^ can be eliminated adiabatically, and the phase 
variable description involving is complete. 

In this paper, we have restricted ourselves to a one-dimensional layered structure. However, the theory is easily 
extended to those in higher dimensions such as a hexagonal structure of disk-shaped domains and a face centered 
cubic structure of spherical domains. Evaluation of nonlinear terms is also straightforward although the calculations 
are more involved. 
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APPENDIX 



Here we summarize some of the steps in the derivation of the coupled set of interface equations of motion in section 



3. From eqs. ( J2 2| ) and (24), the value of v at an interface is given explicitly by 

■jTTyl / d y'Yl[ dx ' 2 i 2 ex P(^( 6 » - x') + iq v (y - y')) 

\ Z1T > J ™ JbZ, Q -t- K 



dx' 3— 2 exp(tq x (bt -x')+iq y {y-y'))] (92) 
q z + n z 



Note that contains the deformations Ctiv) as m ec l- (H|)- 

The zero-th order is obtained by ignoring From (|9^) it is given by 
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(0)± 



a 



/dq f i-rru+w ^ 

?2^)2 / d v'^I} I dx ' i _i_ K 2 ex Pfe(?^ ±w)- iq x x' + iq y (y - y')) 

dx' ^——^ exp(iq x (n£ ±w) - iq x x' + iq v (y - y'))] 

mi+w q + K 

f _^xy^r f dx'^r— — -exp(iq x (n£±w) -iq x x') 

(m+l)£— u; ^ 

dx' — 2 exp(ig x (n.£ ± — icfoa; )] (93) 

mt+w Qx ' K 

The integral over x' in (^) is readily carried out. By using the Poisson summation formula (p9|), one can obtain the 
zero-th solution (|25|). 

In order to calculate the first order correction Vj , one needs to expand and ^ contained in 6^ and 6^, 
respectively, in ( |92| ) . Hence the correction consists of two parts 

prnl-\-w 



/dq pmt-\-w ^ 

— dx' - x exp(iq x (n£±w) - iq x x') 
Z7r m Jml-w Qx + K 



iq x 

ml+w Qx K " 



exp(iq x (n£ ±w) - iq x x')] 



/dcjf f 1 
/ X! 2 T K 2 exp(ig x (nl ± w - ml) + ^(2/ - y')) 

x lC(y')e-^ w - C m (y')e^ w - C n+ i{y')e- iqAl - w) + C+(y')e- iq * w ] (94) 



Evaluation of the term which contains the factor Cn (iO ^ s straightforward, which gives us the first term in (g6|) and 
@. The integrals over q y and y' in the remaining term of ( p^ ) can be easily carried out by introducing the Fourier 
transform (p8|). Using the formula (^), one finally obtains the second term in (|26| ) and (p7|). 
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FIG. 1. Gently deformed layered structure. The dotted lines indicate an equilibrium configuration with w — £/2. 
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